Universality class of triad dynamics on a triangular lattice 
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We consider triad dynamics as it was recently considered by Antal et al. [T. Antal, P. L. Krapivsky, 
and S. Redner, Phys. Rev. E 72 , 036121 (2005)] as an approach to social balance. Here we 
generalize the topology from all-to-all to a regular one of a two-dimensional triangular lattice. The 
driving force in this dynamics is the reduction of frustrated triads in order to reach a balanced state. 
The dynamics is parameterized by a so-called propensity parameter p that determines the tendency 
of negative links to become positive. As a function of p we find a phase transition between different 
kind of absorbing states. The phases differ by the existence of an infinitely connected (percolated) 
cluster of negative links that forms whenever p < p c . Moreover, for p < p c , the time to reach the 
absorbing state grows powerlike with the system size L, while it increases logarithmically with L for 
p > p c . From a finite-size scaling analysis we numerically determine the critical exponents f3 and 
v together with 7, r, a. The exponents satisfy the hyperscaling relations. We also determine the 
fractal dimension df that fulfills a hyperscaling relation as well. The transition of triad dynamics 
between different absorbing states belongs to a universality class with new critical exponents. We 
generalize the triad dynamics to four-cycle dynamics on a square lattice. In this case, again there is 
a transition between different absorbing states, going along with the formation of an infinite cluster 
of negative links, but the usual scaling and hyperscaling relations are violated. 

PACS numbers: 89.75.Fb, 05.40.-a, 64.60.Ak 



I. INTRODUCTION 



Recently Antal et al. Q proposed a so-called triad dy- 
namics to model the approach of social balance 0, OH • 
An essential ingredient in the algorithm is the reduction 
of frustration in the following sense. We assign a value 
of +1 or —1 to a link (or bond) in the all-to-all topology 
if it connects two individuals who are friends or enemies, 
respectively. We call the sign ±1 of a link its spin. If the 
product of links along the boundary of a triad is nega- 
tive, the triad is called frustrated (or imbalanced), other- 
wise it is called balanced (or unfrustrated) . The state of 
the network is called balanced if all triads are balanced. 
The algorithm depends on a parameter p G [0, 1], called 
propensity parameter. It determines the tendency of the 
system to reduce frustration via flipping a negative link 
to a positive one. For an all-to-all topology Antal et al. 
predict a transition from an imbalanced non-absorbing 
stationary state for p < 1/2 to a balanced absorbing 
state for p > 1/2. Here the dynamics is motivated by 
social applications so that the notion of frustration from 
physics goes along with frustration in the psychological 
sense. The mathematical criterion for checking the sta- 
tus of frustration is the same. 

In a recent paper [I| we generalized the triad dynamics 
in two aspects. The first generalization refers to a k-cycle 
dynamics which contains the triad dynamics for k = 3. 
Here it turned out that the main difference comes from 
the difference of whether k is even or odd, since the phase 
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structure is symmetric about p = 1/2 for k even. Even in 
the infinite- volume limit there are only absorbing states, 
apart from the transition point at p = 1/2. The sec- 
ond generalization concerned the network topology from 
all-to-all connections to a diluted network. We studied 
the phase structure as a function of the propensity p and 
the dilution. As it turned out, the diluted k-cycle dy- 
namics can be mapped on a certain satisfiability problem 
in computer science, the so-called fc-XOR-SAT problem 
and socially balanced states in one problem corre- 
spond to all logical constraints satisfied in the fc-XOR- 
SAT problem. In both models we have phases of imbal- 
anced non-absorbing states in the infinite- volume limit, 
separated by a phase transition from phases of balanced 
or absorbing states. In a finite volume one only observes 
balanced states, but as a remnant of the infinite-volume 
phase structure, the time to reach the absorbing states 
differs in a characteristic way. 

In this paper we study triad dynamics on a two- 
dimensional triangular lattice and 4-cycle dynamics 
(called tetrad dynamics) on a square lattice. From the 
interpretation as an approach to social balance, triad or 
tetrad dynamics on a regular topology are not more real- 
istic than on an all-to-all topology. Still, triad dynamics 
shows interesting features in terms of a percolation tran- 
sition if we compare snapshots of frozen states for dif- 
ferent values of p. Also for tetrad dynamics we observe 
a transition between different absorbing states, but the 
description in terms of a percolation transition fails. As 
we shall see, due to the restrictive topology, imbalanced 
triads and tetrads allow only two elementary processes: 
either they diffuse or they annihilate each other. As a re- 
sult, the system always approaches a balanced absorbing 
state; nevertheless we observe a transition as a function 
of the propensity parameter p, this time between dif- 
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ferent absorbing states. The difference is characterized 
by the presence or absence of an infinite cluster of con- 
nected unfriendly (negative) links and for triad dynamics 
by the time to reach these frozen states. The parameter 
p should not be confused with the occupation probability 
of a single bond with a positive sign as it usually used 
in connection with bond percolation. Therefore, here for 
p < p c we observe an infinite cluster of negative links, 
while such a cluster is absent for p > p c . For triad dy- 
namics the time to reach the frozen state grows logarith- 
mically with the system size for p > p c and as a power 
of the size for p < p c . In contrast, for tetrad dynamics it 
grows as a power of the system size only at the transition 
point, while it grows logarithmically above and below p c . 
For triad dynamics we numerically determine the critical 
exponents from a finite-size scaling analysis of the frozen 
patterns. We analyze the probability for a link to 
belong to the infinite cluster that serves as an order pa- 
rameter as well as the order parameter susceptibility and 
the distribution of clusters of finite size s. The critical 
exponents (3, v, 7, a, r fulfill the usual hyperscaling rela- 
tions as well as the fractal dimension df. The critical ex- 
ponents turn out to be new and completely different from 
those of standard percolation in two dimensions [3, |8| ■ 
The outline of the paper is the following. In sectionlTTlwe 
define the triad dynamics on a triangular lattice. In sec- 
tion II I II we analyze the finite- volume dependence of the 
frozen states as a function of the propensity parameter p 
and present the results for the critical exponents as well 
as for the fractal dimension. Section llVl deals with tetrad 
dynamics on a square lattice for which the hyperscaling 
relations are violated. In section [V] we summarize the 
conclusions. 



II. TRIAD DYNAMICS 

Our dynamical system is defined on an undirected 
graph (network) composed of nodes and links. Each 
link between the nodes i and j takes spin values 

a (i-j) = ~1 or = +1 if the nodes i and j are "en- 

emies" or "friends", respectively. A triad [A] (i,j,k) is 
characterized by the values assigned to its three links 
(i,j), (j, k) and (k,i). We have four types of triads, de- 
pending on the number of negative links they contain in 
their boundary: Ao, Ai, A2 and A3 , where the subscript 
stands for the number of negative (or "unfriendly") links. 
We use the standard notion of social balance as proposed 
in 0, Q and apply this notion to triads. The sign of a 
triad is defined as the product of the spins assigned to the 
links of the triad. A triad is considered as "balanced" or 
"unfrustrated" if its sign is positive, otherwise it is called 
"imbalanced" or "frustrated" . The triads Ao (all friends) 
and A 2 (two friends have the same enemy) are balanced, 
while the triads Ai and A3 are imbalanced. The network 
itself is called balanced if and only if all triads belonging 
to the network are balanced. 

As in Q, we perform a local unconstrained dynamics 
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Figure 1: Local dynamics of one update event. The imbal- 
anced triads are represented as filled triangles, while the bal- 
anced ones as empty triangles. The shared link is represented 
as bold, this is the link involved in the update event. A) When 
the update event flips the link shared by one imbalanced and 
one balanced triad we have diffusion. B) When the update 
event flips the link shared by two imbalanced triads we have 
annihilation. 

in order to reduce the frustration of the network. As it 
turns out, the local algorithm always drives the network 
to a fully balanced state without frustrated triads, but 
the time it needs for reaching the frozen state depends on 
the choice of parameters. At each update event one triad 
is selected at random. If the selected triad is balanced 
(type A or A 2 ) nothing happens. If the selected triad is 
imbalanced (type Ai or A3) it is updated into a balanced 
one by flipping one of its link. In particular a triad Ai is 
changed with probability p into a triad Ao (flipping the 
only negative link of the triad Ai), and it is changed with 
probability 1 — p into a triad A2 (choosing at random one 
of the two positive links belonging to the triad Ai and 
inverting the spin to a negative value). A triad A3 is 
changed into a triad A2 with probability 1, choosing at 
random one of its three negative links and reversing its 
spin to a positive sign. We summarize the updating rules 
of the local algorithm in the following scheme: 

A «- Ax > A 2 <- A 3 . (1) 

p 1— p 1 

One time unit has passed when the number of single up- 
date events equals the total number of links M of the 
network. 

We study the triad dynamics on two-dimensional trian- 
gular lattices with periodic boundary conditions. We 
characterize a triangular lattice using its linear size L. 
The total number of sites in the lattice is N = L(L — 1). 
The term L — l results from the periodic boundary condi- 
tions. The total number of links of the lattice is M = 3N, 
while the total number of triads is = 2N. In par- 
ticular each link is shared only by two nearest-neighbor 
triads, so that the triad dynamics cannot increase the to- 
tal number of imbalanced triads: a single update changes 
one selected imbalanced triad into a balanced one, while 
it modifies the other triad, sharing the same updated 
link, either from balanced to imbalanced [Figure ID^-)] or 
from imbalanced to balanced [Figure ^p)]. The former 
we call diffusion, because the imbalanced triad diffuses, 
the latter annihilation between two imbalanced triads. 
Obviously, for a finite-size system we always observe a 
frozen configuration as the stationary state (see Figure 
|2J), independently on the initial configuration. From now 
on we focus on these frozen configurations and study their 
geometrical properties by using the standard tools of per- 
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eolation theory |7j. It should be noticed, however, that 
the dynamical parameter p of the triad dynamics is very 
different from the occupation probability as it is defined 
in percolation theory, where it is also called p. 



III. NUMERICAL SIMULATIONS TO 
DETERMINE THE CRITICAL THRESHOLD 
AND THE CRITICAL EXPONENTS 

In order to study the geometry of the frozen configura- 
tions, we consider the distribution of negative links on the 
lattice. In particular we numerically compute the prob- 
ability that a negative link belongs to an infinite cluster 
Poo = M^/M, as the ratio of the number of links be- 
longing to the largest cluster of connected negative links 
M~ to the total number of links M. 




Figure 2: Typical frozen configurations for a triangular lattice 
with periodic boundary conditions, with linear size L — 129 
and for different values of p: A) p — 0.44 , B) p — p c = 0.4625 
, C) p = 0.48. In the plots only negative links are shown. 



Critical propensity parameter p c 

Let us first determine the critical point p c . Binder's 
cumulant , defined as the 4-th order reduced cumulant 
of the probability distribution 



U L = 1- 



(P 4 



(2) 



allows us to determine the critical point p c without tun- 
ing any parameter. In Eq.J2J) (•} stands for the average 
over all realizations of the distribution. It is known that 
the Binder cumulant should satisfy the scaling relation 



9] 



U L = U (p- Pc )L 



l/u 



(3) 



with £/(•) a universal function. From Figure [3] it is 
obvious that the Binder cumulants for different linear 
sizes of the lattice L have a common intersection at 
p = 0.4625(5), so that we use p c = 0.4625 as the critical 
propensity parameter in the following analysis. The nu- 
merical results of Figure [3] are extracted from numerical 
simulations for L = 17, 33, 49, 65, 97, 129, 193 and 257. 
The averages are taken over several frozen configurations 
as they are reached as absorbing states of the triad dy- 
namics, starting from initial conditions where each link 
is randomly assigned a value of +1 or —1 with the same 
probability 1/2. The number of realizations is 10 4 for 
values of L up to 65 and 10 3 for larger values of L. 



Critical exponent v 

In order to extract the value of the critical exponent 
v that characterizes the divergence of the correlation 
length £ in the vicinity of the critical point, we con- 
sider the absolute value of the first derivative of the 
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Figure 3: Binder cumulant Ul as a function of the propensity 
parameter p. The main plot shows Ul in the vicinity of the 
critical point. The numerical results for different linear lat- 
tice sizes L have a common intersection at p c ~ 0.463. The 
inset shows the log-log plot of \dV 'l / 'dp\ p . From the relation 

\dU L /dp\ pc ~ L 1/V we find l/v = 0.64(7) (solid line). 



Binder cumulant calculated at p Cl \dU L / dp\ p ^. We ex- 
pect from Eq.Q to have \dU L /dp\ pc ~ L 1 ^. This re- 
lation is actually satisfied for l/v = 0.64(7) as it is 
seen in the inset of Figure |3 from which v — 1.6(2). 
For the numerical estimate of the derivative of Ul we 



use \dU L /dp\ 



E?=i q a L, q p, 



9-1 



where the coef- 



ficients a Ll q are obtained by interpolating the cumulant 



Ul with a polynomial Ul 



<lL,q P H 



Here we 



choose Q = 5 and restrict the interpolation to the inter- 
val [0.46,0.465]. 



Critical exponent /3 

According to percolation theory, the probability P a 
satisfies the following finite-size scaling relation 



(4) 



Figure 4: Probability P x that a bond belongs to the largest 
cluster as a function of p and for different linear sizes of 
the lattice. We perform a finite-size scaling for determin- 
ing the critical point p c — 0.4625(5) [see the zoom around 
the critical point in the inset A)] with the critical exponents 
fifv = 0.297(3) (solid line) as is shown in the inset B), where 
Poc at the critical point is plotted as a function of the linear 
size of the lattice L. The numerical simulations are the same 
as in Figure |21 




Figure 5: Same data as in Figure 2] but now the abscissa is 
rescaled according to (p — p^L 1 ^ with l/v — 0.64. The inset 
shows a zoom around zero. The numerical simulations are the 
same as in Figure |3 



where p c is the critical value of the propensity param- 
eter p for which the phase transition occurs. (In the 
infinite- volume limit we have = 1 for p < p c , while 
Poo = for p > p c .) P(-) is a universal function. A plot 
as a function of p and for different values of L 
is shown in Figure For a value of f3/v — 0.297(3) [see 
inset B)] all curves have an intersection inp c = 0.4625(5) 
as it is more obvious from the inset A). Using for v the 
value calculated so far, we obtain (3 — 0.46(5). 
If we rescale the abscissa as (p—p^L 1 ^ with p c — 0.4625 
and \ jv = 0.64, all curves for different values of L col- 
lapse into one (see FigurelSJ). This is true especially close 
to zero as the inset of Figure [S] clearly shows. 



Critical exponent 7 

Moreover, in Figure we plot the susceptibility 

X = M[{Pl)-{P oa f] , (5) 

in which we only show a zoom around p — p c , while in 
the inset we plot the value of %, at the critical point, as 
a function of the linear size of the lattice L. From the 
inset we find 7/1/ = 1.28(3), because \ ~ -Zj 7 ^ at p c , 
therefore 7 = 2.0(3). The susceptibility should satisfy 
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Figure 6: Susceptibility \ of as a function of p. In the 
main plot we zoom the region around the critical point p c , 
while in the inset we plot x a ^ the critical point as a function 
of L leading to 7/f = 1.28(3) (solid line). The numerical 
simulations are the same as in Figure 
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Figure 7: finite-size scaling for the susceptibility \- The data 
are the same as in Figure|S| In the inset we zoom the function 
\L around the critical point p c , while in the main plot we 
rescale the abscissa as (p — p c )L}^ v . The critical exponents 
are chosen as i/v = 1.28 and \jv = 0.64. 



the finite-size scaling relation 

X = L^xUp~ Pc )L 1 /» 



(6) 



where again x(-) is a universal function. This relation is 
perfectly satisfied, cf. Figure [7| 



Critical exponents r and a 



Figure 8: The main graph shows the distribution of the cluster 
size n a (p c ) at the critical point p c and for L = 257. The 
distribution is extracted from 10 3 frozen configurations. As 
expected this distribution follows a power law n s ~ s~ T , the 
solid line plotted here corresponds to r = 2.19. The inset 
shows the second moment of the distribution of the cluster 
size M2 at the critical point as a function of the linear size of 
the lattice. M2 increases as a power of the linear size L of the 
lattice with exponent y/v — 1.28(2) (solid line). 



theory, n s should satisfy 



= s T h [(p - p c 



(7) 



where r and a are critical exponents and n(-) is a uni- 
versal function. We can determine the critical exponent 
r (also called Fisher exponent) by plotting the distribu- 
tion of the cluster sizes s at the critical point p c . This 
distribution is calculated in Figure for L = 257. The 
distribution fits with a power law s~ T , here plotted as 
solid line with r = 2.19. We checked that the same ex- 
ponent fits also for smaller values of L. 
In Figure El we numerically determine the critical ex- 
ponent a. The figure shows the plot of the ratio 
n s (p)/n s {p c ). Using the rescaled variable (p — p c )s a as 
abscissa, all the curves corresponding to different val- 
ues of the cluster size s ( 2 4 < s < 2 5 , 2 5 < s < 2 6 , 



2 b < s < 2\ etc. 



collapse to a single function for 



a = 0.41, as expected, see e.g. 



Check of the critical exponent 7 

From the cluster distribution it is possible to verify 
the value of the critical exponent 7 that was calculated 
before. It is known that the second moment of the dis- 
tribution of the cluster size 



M 2 = J2 s2n * 



(8) 



Furthermore, we consider the probability distribution 
of having n s clusters with s negative links. n s is given 
by the ratio of the number of clusters of size s to the 
total number of clusters. As it is known from percolation 



should scale according to Mi ~ L 1 ^ at the critical point 
p c . The numerical value found for the ratio 7/V = 1.28(2) 
is consistent with the former one obtained via the suscep- 
tibility (see the inset of Figure [SJ. 
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Figure 9: Numerical test of the validity of equation J7J. We 
plot the ratio n B (p)/n s (p c ), where n 3 is the distribution of the 
cluster size, as a function of the rescaled quantity (p — p c )s" 
with a chosen as 0.41 from HIIII knowing j3 and 7. The main 
plot refers to L — 257, the inset to L = 129. 



Fractal dimension df 

The percolating cluster can be further characterized in 
terms of a fractal dimension df. The fractal dimension df 
is easily computed using a box counting method We 
count how many negative links M(r) belong to a trian- 
gular box with r links per side. The measurement is per- 
formed for L = 257 at the critical point p c . We analyze 
10 3 frozen configurations, for each of them we consider 
10 2 different boxes. The mass M(r) per box is plotted 
in Figure El As expected we observe the crossover phe- 
nomenon 



M(r) 



r d f 



, if r < £ 
, if r > £ 



(9) 



where d = 2 in this case. The numerical results fit with 
a power law with df — 1.703 and d — 2, respectively. 



Time to reach the frozen states 

So far we characterized the geometrical properties of 
the final absorbing configuration. As a next step we fo- 
cus on the dynamical features and determine the time 
the system needs to reach the final frozen configurations. 
In Figure IllA) we plot the ratio of the variance ctt and 
its average value T as a function of p. It is interesting 
to notice that ctt/T has a maximum around p c . Below 
the critical point p c , the relaxation time T is governed 
by a power-law relation in the linear size of the lattice L: 
T ~ L z . Our numerical analysis reveals z — 2.36(1) at 
p c and z = 2.24(1) at p = 1/3 [see Figure Hlk)]. Further- 
more, above the critical point p c , T grows logarithmically 
with L: T ~ log (L) [see FigurelrTDh where we show the 
L dependence of T for p = 3/4 and for p = 1]. The 



Figure 10: Mass of negative links M inside a triangular box 
with r links per side. The mass is measured on a triangular 
lattice with L = 257 and periodic boundary conditions. Each 
point is given by the average over 10 3 frozen configurations, 
while the number of boxes considered for each r is 10 2 . The 
numerical results fit a power law with df = 1.703 (solid line) 
and d = 2 (dotted line) . 



different size dependence of the time the system needs 
to reach a frozen configuration in both phases can be 
understood in qualitative terms. For p > p c there is a 
lower probability for having negative links. The stable 
configurations forming out of these links are mainly local 
objects like star-like configurations which do not perco- 
late through the lattice [see Figure HP)]. Therefore the 
time to reach such a state characterized by local objects 
depends weakly on the system size, that is logarithmi- 
cally. For p < p c there is a high probability of having 
negative links. However, here large loops that finally 
percolate through the lattice make up the stable config- 
urations [see Figure IflV) and B)]. So the time to reach 
this configuration is more sensitive to the system size, it 
grows like a power of the size. 



Universality class of triad dynamics in two 
dimensions 



We list the critical exponents for the phase transition 
between different absorbing states in Tabled The critical 
exponents satisfy the known hyperscaling relations. For 
example, it is known from percolation theory that the 
critical exponents 0, 7, r and a are related by 



T= 2 



1-0 



and 



7 + /3 



(10) 



(11) 



Calculating r and a from the former equations, using 
the numerically obtained values of f3 and 7, we find r — 
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Figure 11: Average time T needed for reaching a frozen con- 
figuration. The datasets are the same as in Figure |I] In A) 
we plot the standard variance ctt over the average value T as 
a function of p and for different values of L. The peaks of 
this ratio are all around the critical point p c . In B) we show a 
log- log plot of T as a function of the linear size of the lattice L 
at the critical point p c (full circles) and at p = 1/3 (open cir- 
cles). It implies that T ~ L z for p < p c . We have z = 2.36(1) 
(solid line) at p c , while z = 2.24(2) at p = 1/3 . In C) the 
plot of T versus the logarithm of L for p — 1 (full circles) and 
for p — 3/4 (open circles) shows that T ~ log (L) for p > p c . 
The numerical simulations are the same as in Figure |H] 



Pc 


0.4625(5) 


p 


0.46(5) 


7 


2.0(3) 


V 


1.6(2) 


a 


0.41(6) 


T 


2.19(1) 


df 


1.703(3) 


z 


2.36(1) 



Table I: Critical threshold and critical exponents numerically 
determined for the absorbing configurations of triad dynamics 
on two-dimensional triangular lattices. 



2.19(1) and a = 0.41(6). Both values are consistent with 
Figures Eland El 

The fractal dimension is related to the ratio (3/v by the 
hyperscaling relation 

d f = d--. (12) 
v 

From this relation we find df = 1.703(3) in agreement 
with Figure ITUl 

The errors of p c and of the critical exponents as indi- 
cated in Table |U arise as follows. The critical point p c is 
directly estimated from Figure |3 verified in Figure Q] and 
supported by Figure ^] The only source of error here 
is given by the step size used for varying the propensity 
parameter p. Of course this step size is related to the 
parameters of the simulations. It can be decreased by 



increasing the system size and the total number of simu- 
lations. 

The critical exponent (3/v of Figure0j3) , \ jv of Figured! 
7/1/ of Figures0and|H| and z of FigureHTlhave errors due 
to the linear fit in a double-logarithmic plane. In princi- 
ple one should also account for the propagation of the er- 
ror entering the value of the critical point p c — 0.462(5), 
which we here have neglected. 

In order to check the validity of the hyperscaling relations 
we evaluate the corrections to all derived quantities (t, a 
and df) by using the standard formula for error propa- 
gation. The hyperscaling relations are then said to be 
satisfied if they hold within these error bars. 
To our knowledge the critical exponents calculated so far 
are new [?l l§l flfj|. Therefore, the percolation transition 
between different absorbing states of triad dynamics on 
two-dimensional lattices can be described by standard 
percolation theory, but the transition seems to belong to 
a new universality class. 

IV. TETRAD DYNAMICS 

The notion of social balance can be extended to any 
geometric figure with k links (in graph theory denoted 
as fc-cycles) 0, H|. As it turned out in 5] for an all- 
to-all topology and a diluted topology, generalizing the 
triad dynamics to a fc-cycle dynamics leads to qualitative 
differences in the phase structure of which the main dif- 
ferences are due to k being even rather than being larger 
than three. In this section we therefore focus on the case 
of k = 4 and call the 4-cycles tetrads defined on square 
lattices. According to the number of negative links be- 
longing to a particular tetrad, we distinguish five types 
of tetrads: Do , Di , D2 , D3 , LI4 . The updating rules of 
the unconstrained dynamics for triads of Q are naturally 
extended to tetrads in the following way: 

□0 <- Di ► D 2 «- D 3 ► D 4 . (13) 

p l— p p i—p 

The local tetrad dynamics is then applied to square lat- 
tices with periodic boundary conditions. L denotes the 
linear size of the lattice, i.e. the number of sites per row 
or per column. The total number of sites of the lattice 
is N — L 2 , the total number of links M = 2N, and the 
total number of tetrads N[j — N. 



Similarly to the case of triad dynamics on two- 
dimensional triangular lattices, at a first glance there 
seems to be a percolation transition also in case of tetrad 
dynamics. The critical point p c should be equal to 
Pc = 0.5 due to the symmetry of the system (|13|) un- 
der the simultaneous transformation cr, — » — <7j Vi and 
p — » 1 — p (see Figure 1121 . 

However, differently from the triad dynamics, no scal- 
ing relations hold in the vicinity of the transition at 
which again an infinite cluster of negative links forms. 
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Figure 12: Typical frozen configurations for a square lattice 
with periodic boundary conditions, with linear size L — 128 
and for different values of p: A) p = 0.48 , B) p — p c — 0.5 , 
C) p — 0.52. In these plot are shown only the negative links. 
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Figure 13: Probability Poo that a link belongs to the largest 
cluster as a function of p and for different linear sizes of the 
lattice. The ratio of critical exponents j3/v = 1.05(2) at the 
critical point p c — 0.5 as is shown in the inset. 




Figure 14: Mass of negative links M inside a quadratic box 
with r links per side. This measure is used on a square lattice 
with L — 256 and periodic boundary conditions. Each point 
is given by the average over 10 3 frozen configurations, while 
the number of boxes considered at each r is 10 3 . The solid 
line has the slope 0.95, a value that is expected by the former 
measure of /3/v. Actually it is not consistent with the direct 
fit of the data from which we find df — 1.76(2) (dashed line). 
The dotted line has slope d — 2. 



The hyperscaling relations are violated. For example, 
the relation Ijl2(l does not hold because one can mea- 
sure (3/v = 1.05(2) (Figure I13|) and expect to have 
df = 0.95(2), while it is directly seen from Figure ITU 
that d f = 1.76(2). 

Tetrad dynamics resembles triad dynamics. Each un- 
stable tetrad either diffuses, or two unstable tetrads an- 
nihilate if they meet. Finite-size systems always reach a 
frozen configuration within a finite time, so we can mea- 
sure this time T. As we can see from Figure ITSR ) . we 
find a symmetric behavior around the critical point p c . 
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Figure 15: Average time T needed for reaching a frozen con- 
figuration. The simulations are the same as in Figure FTTH I n 

A) we plot the standard variance <tt over the average value 
T as a function of p and for different values of L. The peaks 
of this ratio are all in the vicinity of the critical point p c . In 

B) we show a log-log plot of T as a function of the linear size 
of the lattice L at the critical point leading to T ~ L z with 
z = 2.35(1). In C) the plot of T versus the logarithm of L for 
p = 1 shows that T ~ log (L). 



At the critical point T scales according to T ~ L z with 
z = 2.35(1) [see Figure IPoE?)]. this value is actually con- 
sistent with the one found for the triad dynamics. Away 
from the critical point T ~ log(L) [see Figure USD)], 



The driving force in triad dynamics is the reduction 
of the number of frustrated triads. A state of zero frus- 
tration is called a state of social balance. Imposed on a 
triangular lattice, frustrated triads can diffuse or annihi- 
late each other. Depending on the value of the propen- 
sity parameter the final absorbing state can be charac- 
terized by the absence (p > p c ) or presence (p < p c ) of 
an infinite cluster of negative links. The time to reach 
the frozen configurations scales with the system size in a 
way that further characterizes the phases: it scales log- 
arithmically for p > p c and power-like for p < p c . The 
critical exponents v, [3, 7, r, and a as well as the frac- 
tal dimension d f satisfy hyperscaling relations within the 
error bars. The values of these exponents seem to char- 
acterize a new universality class. The essential difference 
that we observe for tetrad dynamics on a square lattice 
is the symmetry between the absorbing states in the dif- 
ferent phases. Again, an infinite cluster of negative links 
forms for p < p c , but the time to reach the frozen config- 
urations shows the same dependence on the system size 
in both phases, with the only exception at the transition 
point. The percolation picture breaks down in the sense 
that no scaling and therefore no hyperscaling relations 
are satisfied. 
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